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Abstract 

Exact integrability of the Sasa Satsuma equation (SSE) in the Liou- 
ville sense is established by showing the existence of an infinite set of 
conservation laws. The explicit form of the conserved quantities in terms 
of the fields are obtained by solving the Riccati equation for the associ- 
ated 3x3 Lax operator. The soliton solutions, in particular, one and two 
soliton solutions, are constructed by the Hirota's bilinear method. The 
one soliton solution is also compared with that found through the inverse 
scattering method. The gauge equivalence of the SSE with a generalized 
Landau Lifshtz equation is established with the explicit construction of 
the new equivalent Lax pair. 
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I. Introduction 

Nonlinear Schrodingcr equation and its various generalized versions (higher or- 
der nonlinear Schrodinger equation) is well known in describing various physical 
phenomena [Q . A common property in all these physical systems is the appear- 
ance of solitons, as a result of a balance between the nonlinear and dispersive 
terms of the wave equations. With the advancement of experimental accuracy, 
solitons having more complicated dynamics can also be detected and observed 
now ||] . The Sasa Satsuma equation || 

*Qt + \Qxx + \Q\ 2 Q + ^{Qxxx + G\Q\ 2 Qx + 3\Q\xQ) = (1) 

describing the evolution of a complex scalar field, is an example of such a system, 
whose soliton solutions have been obtained through inverse scattering method 
(ISM) in §. 

A limited class of soliton bearing equations exhibits further interesting prop- 
erties and belongs to the exclusive club of integrable systems. The most promi- 
nent definition of integrability is the integrability in the Liouville sense, i.e., the 
existence of a set of infinite numbers of conserved quantities in involution Q , 
which can be considered as the action variables. This criterion of integrability 
is extendable also to the quantum case. The Lax pair associated with the model 
is usually a sign of such integrability, while the Painlev e singularity analysis |^| 
is supposed to be a direct test of integrability for the given equation. 

It should be mentioned that the Lax pair for the SSE as well as its one soliton 
through ISM were found in || , while the Painleve analysis for the equation was 
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carried out in ||. However, extraction of the higher conserved quantities for 
the Sasa Satsuma system and thus establishing the intcgrability of the whole 
hierarchy in the Liouville sense remained unexplored. A possible reason of this 
may be the difficulty involved due to the unusual 3x3 matrix form of the Lax 
operator, associated with the SSE. 

Our objective is, therefore, to find the Riccati equaton for the 3x3 Lax op- 
erator, associated with the SSE and consequently, to obtain the whole hierarchy 
of conserved charges in a systematic way. This constructiom will be somewhat 
involved due to the extended form of the Lax operator. 

For investigating SSE from a different view point, we further find the explicit 
soliton solutions of the equation through Hirota's bilinear method. This is a di- 
rect and much more effective method compared to ISM for obtaining the soliton 
solutions, since it does not require the knowledge of the Lax pair. Moreover, 
the construction of the r function becomes straightforward in this method. 

One should recall in this context another interesting fact about the NLS 
equation that it is gauge related to the well known Landau Lifshitz equation 
(LLE) 0. This equivalence can also be established through the space curve 
method B. It is, therefore, natural to ask what is the gauge equivalent equation 
to the SSE. Though such an equivalent system has been discovered by the 
space curve method in H], we complete the investigation for SSE by showing 
the equivalence of it through an explicit gauge tranformation, which not only 
reproduces the generalized LLE (GLLE), but also constructs the associated Lax 
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pair for the GLLE. 

The organization of this paper is as follows. In section II we study the soliton 
solutions of SSE by Hirota's bilinear method. We compute explicitly the one 
and two soliton solutions and compare our result of one soliton solution with the 
known one S. In section III, we construct the related Riccati equation using 
the 3x3 Lax operator of the SSE and subsequently find the infinite number of 
conserved quantities through the recursion relation. The time invariance of the 
conserved quantities is checked directly by using the evolution equation. This 
proves the Liouville integrability of the SSE. The section IV, provides the gauge 
equivalent generalized LLE and gives the associated new Lax operators in the 
explicit form. Section V is the concluding one. 

II. Soliton Solutions Through Hirota's Method 

Let us begin with the SSE which through a change of variable and a 
Galelian transformation: 

Q(X,T) = u(x,t)exp{ie(x+^)} 

6 

T = t (2) 
X = x+U 

may be simplified to the form 

u t + ^(u m + 6 | u | 2 u x + 3(1 u \ 2 ) x u) = (3) 

This is an example of a complex modified KdV type equation and goes to mKdV 
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for the real valued field. The associated spectral problem can be studied through 
the pair of linear equations 

^ x = U(x,t,\)V (4a) 
#, = V(x,t,X)9 (4b) 

where U(x, t, A) and V(x, t, A) are 3x3 matrices and A is the spectral parameter. 
The explicit form of U and V may be given using the result of as 

U = -iXE + A (5a) 
V = -i4eA 3 £ + 4e(A 2 - | u | 2 )A 

-i2eAE(A 2 - A x ) - eA xx + e(A x A - AA X ) (5b) 

with 




/ u \ 
A= u* J 
\ -u* -u / 

Compatibility of (|a) and (|b) leads to SSE ©, which can be shown easily by 
using the following properties of S and A matrices 

{S,A} = 
S 2 = 1 
A 3 + 2|u| 2 A = 

Note that SSE in the form (|^), though suitable for studying inverse scattering 
technique, is not convenient for casting it into Hirota's bilinear form. On the 
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other hand, the higher order nonlinear Schrodinger equation form ([!]) for SSE is 
more suitable for this purpose. Now, in order to write (|l|) in the bilinear form, 
we make the transformation 

Q{X,T) = G{X,T)/F{X,T) (6) 

where, G is complex and F is real. Consequently, in these new variables, we 
have the following set of equations 

(iD T +lD 2 x + l-D 3 x )G.F = (7a) 
D X F.F = AG* .G (7b) 
(l-jD x )G*.G = (7c) 

which follow from (|l|). Dt, Dx, Dxx etc. in (Q) are Hirota derivatives (0) 
belongs to a new class of bilinear equations, whose general form would be of the 
type 

B{D X ,D T ,---)G.F = (8a) 
A(D X ,D T ,---)F.F = C(D X ,D T ,---)G*.G (8b) 
£(l- D X ,D T ,---)G*.G = (8c) 

The additional bilinear equation (^c) involving G*G imposes one further con- 
dition on the complex parameter P (shown below), which is absent in other 
examples of higher order nonlinear Schrodinger equations |h]] . 

For obtaining one soliton solution of SSE d|), we choose G and F in the 
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following form 

G = iexp(?7) (9a) 
F = 1 + K exp(r? + rf) (9b) 

where, L is a complex c-number parameter and 

r] = PX + flT + (10) 

with P, f2 are in general complex parameters. Substituting the expressions ([)]) 
for G and F in ||), we see that G and F are the solutions of (^) provided the 
following relations hold 

1 P 3 

iQ + -P 2 + i— = (11a) 

2 oe 

A' = ^r am 

P-P*=2ie (11c) 
Where the complex parameter, P is of the form 

P = H + ie (12a) 



(11a) is nothing but the dispersion relation and (lib) determines K. In the 
above solution, so far, L is an arbitrary complex parameter. We will see shortly 
that in order to compare our result with the one obtained through the ISM Q , 
the parameter L is to be chosen in a specific form. It follows from the dispersion 



relation (11a) and the expression of P (12a) that fl should be of the form 



2 2 

^ = -MtJ- + f ) - 4 (12b) 

6e 2 3 



Substituting |t2|) and fll2b| ) in @ and using ( |llb| ) the one soliton solution in 



the explicit form becomes 

L exp 77 



Q(X,T) 



1 + K exp (77 + rj* ) 
Xexpj^ + ^X-^ + f +7^} 
l + ^exp{2 Ai X-(f + M e)T} 

To compare (O) with that of ISM, we choose L as 



(13) 



L = 2fiaq>(-ieX {1) - ^X (0) ) (14) 
which reduces (|l3|) to the form 

JLexp{ie(X-fT-X«)} 

Q(X T) = - (15) 

cosh {fiX - + |)T - M* (0) ) + log V2) 

and this is in agreement with the ISM result of ||. 

Two Soliton solution of the Sasa Satsuma equation may be obtained, follow- 
ing Jl(| , by choosing G and F of the form 

G = L(cxprji + Mexp (771 + tja)) (16a) 

F = 1 + A; exp (771 + 77^ + exp {r) A ) + kMM* exp (771 +77^+ 77,4) (16b) 

where, 771 is complex as in the case of one solition, but 77^4 is chosen to be real. 
Substituting F and G (|l6| ) in the bilinear forms (0), we observe that 771 soliton 



satisfies the same dispersion relation as (11a), whereas the dispersion relation 
for 77^4 becomes P\ = 0. The degree 2 term in (j^a) determines the parameter 
M to be unity. Once again following Jl(| , we define the degree of a term by the 
number of 77's present in the exponent. Now the degree 2 and 4 terms in (Qb), 



yield the value of k as in (lib). However, because of the additional bilinear 
form (0c) in our case, we obtain one more relation, e(P* — Pi) = 2i. Note 
that, in general for complex bosonic systems having complex parameters, the 
two soliton solutions may have some relation, which is analogous to the three 
soliton conditions |ic}| . But, the choice of the real parameter t/a, makes the 
three soliton condition trivial in our case. A more general choice of F and G 
for two soliton solutions will give such a nontrivial condition. 

III. Riccati Equation and The Conserved Quantities 

To show the Liouville integrability, i.e. the existence of an infinite number 
of conserved quantities related to SSE (|J), we first write the associated Riccati 
equation. Since the Lax operators in this case are 3x3 matrices, the Riccati 
equation becomes more complicated, though tactable. Let us write the auxiliary 
field ^ in the component form as 

/ Xlx \ 

* = Xix (17) 
\ X-ix J 

Substituting (^) in (Qa), we get a set of three coupled equations 

Xlx = -i\xi+ux3 (18a) 
X2x = -«Ax2 + u 2 \ 3 (18b) 
X3x = i^X3 - u*xi - UX2 (18c) 

Now expressing ( |lg| ) in terms of T\ = (xi/X3) an( i T 2 = (X2/X3)i an( i elimi- 



8 



nating xi, X2 and %3, one obtains 



X 



u — 2iAri + u*T\ + VT1T2 



(19a) 



r 2 



u* - 2i\T 2 + uY\ + u*T 1 T 2 



(19b) 



The first order nonlinear coupled equations (|l9|) for Ti and T 2 are the Riccati 
equations in our case. Notice that neither the integral of I\ nor of r 2 , plays the 
role of generating functions for conserved quantities, but a suitable combina- 
tion of them does. The infinite number of conserved quantities (Hamiltonians), 
i?2n+i; n = 0, 1, 2, can be obtained from ( |l9| ) by identifying 



a(A) = exp(-iXx) \ x ^oo * 3 (oo, A) = exp{- / (u*Ti + uT 2 )dx} (20a) 



We will see that Hamiltonians with odd indices only survive, while the terms 
with even indices become trivial. This property is similar to that of the real 
KdV or the modified KdV equation. 

We may look for series solutions of ( |l9|) by assuming Ti and T 2 in the form 




where H 2n+1 are related to a(A) as 



oo 




(20b) 



OO 




(21a) 



OO 




(21b) 



which yield the following recursion relations from (19i) and (^l| 
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± 2i 



n+l 



2*Cn+2 — ~(^n+l)x + ^ ( u * ^L^n-m+l + "C^n-m+l)- (22) 
m=0 

Similarly (|l^b) and (|2l]) determine another, though quite similar, set of recur- 
sion relations given as 



n+l 

2*C^ +2 = — (Cn+\)x + ^ ("^M-m+l + u * ^m^n-m+l) (23) 
m— 

Inserting the expressions of C\ and C„, thus obtained through the recursion 
relations (p2[),(p^) in (pl|), we get from (20a,b) the explicit form of all conserved 



quantities, H^n+i. These expressions are the integrals taken over the functions 
of the fields u and u* and their derivatives. The first few conserved quantities 
of the infinite set of the SSE hierarchies are given by 



Hi = I u*udx (24) 
H 3 = J (-u* x u x + 2\u \ A )dx (25) 
#5 = [« x u xx - 4 I u | 2 u* x u x + ((I u \ 2 ) x f + 8 | u f)dx (26) 



We have checked explicitly by using equations of motion (g) that Hi, H3 and 
i?5 are, indeed, the constants of motion. 
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IV. Generalised Landau Lifshitz Type Equation as The Gauge 
Equivalent System 

We now show an interesting connection between the SSE in the form (^) 
and the generalized Landau Lifshitz type equation, by exploiting the gauge 
equivalence of the Lax pairs of these two dynamical systems. The procedure is 
similar to that between the NLS and the standard Landau Lifshitz equation Q . 
Under a local gauge tranformation, the Jost function, ^(x,t,X) changes as 

*(x,t,\)=g- 1 (x,t)9(x,t ) \) (27) 

where g(x, i) — ^{x, t, A)|a=0j may be taken as an element of the gauge group. 
Consequently, the Lax equations (J4J) under this gauge transformation (|27|) be- 
come 

# x = U(x,i,A)§ (28a) 
$ t =V(a:,t,A)* (28b) 
where U and V are the new gauge transformed Lax pair, given by 

U(x,t,\) = g-\V-V )g (29a) 
V(x,t,X) = g-\V-V )g (29b) 

with U = U| A=0 = 9x(x,t)g^ 1 (x 1 t) and V = V| A=0 = g t {x,t)g^ 1 (x 1 t). 
We may identify the spin field of the Landau Lifshitz type equation as 

S = g{x,t)- l Vg{x,t) S 2 = l. (30) 
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With this identification, the gauge transformed Lax pair ( p9| ) can be expressed 
in terms of the spin field S ( |30| ) and its derivatives only, yielding 

U = -iXS (31a) 
V = -4ie\ 3 S + 2e\ 2 SS x +ie\(S xx + ^SSl) (31b) 

In deriving ( |3lj ) one has to use the following important identities 
SS X = 2g~ 1 Ag, SS 2 = -Ag~ 1 HA 2 g and S xx + SS 2 = 2g- 1 H\ x g. 
The zero curvature condition of (|3l|); 

U t -V x + [U,V] =0 (32) 

leads to the generalised Landau Lifshitz type equation 

3 

St + £S XXX + —c(S x + SS XX S X + SS X S XX ) = (33) 
with S e SU(3)/U(1). 

V. Conclusion 

In this paper, we have bilinearised higher order nonlinear Schrodinger equa- 
tion, vis a vis the SSE following Hirota's method. The Hirota's method is an 
effective and important method to obtain multi soliton solutions. We have found 
explicitly one and two soliton solutions and recovered the ISM result of one soli- 
ton solution from that of Hirota's method after a specific choice of the parameter 
involved. The result related to higher soliton solutions are more complicated 
and will be given elsewhere. It is found that the SSE falls under a new class of 
bilinear forms. 
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The linear problem of SSE is a nontrivial one in the sense that the Lax 
operator corresponding to this dynamical equation is a 3 x 3 matrix, which 
makes the related Riccati equation more involved. However, by solving such 
coupled Riccati equations we are able to compute explicitly the infinite number 
of conserved quantities through the recursion relations and to show that terms 
with odd indices only contribute to the conserved charges like KdV and mKdV 
systems. This result establishes explicitly the integrability of the SSE in the 
Liouville sense. Finding out the Poisson bracket structures among the dynamical 
fields and consequently revealing the explicit form of the hierarchy of SSE from 
the conserved charges obtained here will be an interesting future problem. 

The gauge equivalence of the GLLE with the SSE has been established here. 
This equivalence not only gives a direct relationship between the fields, which 
would help to find the soliton solutions of the GLLE using those of the SSE, 
but also yields explicit Lax operators for GLLE from which one would be able 
to extract the related higher conserved quantities in the similar way following 
the present results of SSE. 
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